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A  COMPUTER  SIMULATION  MODEL  OF  OUTDOOR 


RADIANT  LIGHTING 


INTRODUCTION 


The  following  simplified  problem  is  considered.  An  observer  standing 
on  a  flat,  infinitely-extending  ground  surface  is  sighting  with  a  telescope 
in  the  direction  of  a  panel-type  target.  The  ground  layer  is  assumed  to  be 
covered  with  a  material  of  uniform  composition  having  a  diffusive  reflec¬ 
tance,  rfl,  which  is  spectrally  dependent.  The  target  panel  is  a  small  ver¬ 
tical  surface  coated  with  a  diffusively-reflecting  material  having  a  spec¬ 
trally  dependent  reflectance,  rt .  The  panel  is  standing  against  a  sloped 
incline  as  background.  The  background  material  has  a  diffusive,  spectrally 
dependent  reflectance,  rb.  The  exterior  surface  of  the  objective  lens  piece 
in  the  telescope  is  partially  coated  with  a  material  having  spectrally  de¬ 
pendent  coefficients,  t^  and  s^,  for  light  transmittance  and  scattering, 
respectively. 


The  radiant  energy  reaching  the  telescope  from  the  target  and  its 
background  is  a  function  of  the  following  variables: 

1.  The  irradiance  from  the  sky  and  surrounding  ground,  which  im¬ 
pinges  upon  the  target  and  the  background  surfaces. 

2.  The  reflectance  values  of  the  target  and  background  surfaces. 

3.  The  radiance  from  the  horizontal  sky  behind  the  target  and 
background . 

4.  The  atmospheric  extinction  (attenuation)  coefficient  at  the 
ground  level. 


5.  The  observer-to-target  and  observer-to-background  distances. 


The  radiant  energy  reaching  the  telescope  is  attenuated  and  scattered  by 
the  coating  on  the  exterior  surface  of  the  objective  lens.  The  distribution 
of  this  energy  in  the  telescopic  image,  and  the  resulting  luminance  and 
chromatic  contrasts  between  the  target  and  background,  determine  the  visi¬ 
bility  of  the  target. 

The  objective  of  this  report  is  to  establish  equations  for  the  appar¬ 
ent  radiance  at  the  telescope  from  the  target  and  background  as  functions 
of  the  sky  irradiance  and  sun  angles,  and  the  other  physical  and  geometri¬ 
cal  variables.  The  following  sections  establish  sky  radiance  equations  as  a 
function  of  sky  irradiance  and  sun  angles  for  a  clear  sky.  The  ground  ra¬ 
diance,  target  radiance,  background  radiance,  telescope  irradiance  and  the 
apparent  target  and  background  radiances  at  the  telescope  are  next  derived. 
Finally,  the  observer-target  coordinate  system  is  related  to  the  geographi¬ 
cal  coordinate  system  of  the  global  world. 
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DISCUSSION 


Sky  Radiance  Equations 

A  review  of  the  literature  shows  that  extensive  analysis  and  measure- 
aents  have  been  completed  In  the  field  of  sky  lighting.  Chandraskhar  (2) 
used  numerical  techniques  to  develop  tables  listing  the  intensity  and  po¬ 
larization  of  a  pure  air,  planer  atmosphere.  Earlier  work  assuming  vertical 
light  fluxes,  as  well  as  natural  light,  had  been  developed  for  pure  air 
(29)  and  for  overcast  skies  (21).  The  effects  of  the  earth's  curvature  and 
shadow,  important  during  twilight,  have  also  been  considered  (27).  In  ad¬ 
dition,  numerous  measurements  have  been  recorded  of  the  daytime,  twilight 
and  nighttime  skies  (12,  16,  17,  18,  19,  20,  25  and  29). 

We  shall  follow  the  analysis  developed  by  Tousey  and  Hulburt  (29)  to 
derive  equations  for  the  sky  radiance.  They  made  several  simplifying  as¬ 
sumptions:  (1)  the  sky  light  is  natural  and  Stokes'  parameters  for  polar¬ 
ized  light  may  be  ignored,  (2)  the  atmosphere  has  a  uniform  composition, 
(3)  the  atmosphere  is  bounded  top  and  bottom  by  infinitely-extending  paral¬ 
lel  planes,  and  (A)  the  diffused  light  in  the  atmosphere  is  described  by 
opposing  flux  flows  in  the  vertical  direction.  Chandraskhar  (2)  has  shown 
that  the  assumption  of  vertical  light  fluxes  do  not  properly  describe  the 
angular  distribution  of  diffused  light  in  planar  atmospheres.  Furthermore, 
the  assumption  of  planar  atmospheres  is  not  applicable  to  twilight  since 
the  effects  of  the  earth's  curvature  and  shadow  upon  atmospheric  lighting 
is  ignored.  Tousey  and  Hulburt  used  their  analysis  to  develop  brightness 
and  polarization  distributions  for  the  daytime  sky.  Their  results  did  not 
correctly  specify  the  polarization  in  all  sectors  of  the  sky,  but  the 
brightness  values  were  only  in  slight  error  (29). 

Tousey  and  Hulburt  considered  a  planar  atmosphere  of  uniform  composi¬ 
tion  and  constant  thickness,  zt .  The  top  of  the  atmosphere  is  irradiated 
with  sunlight  of  intensity,  Iqx»  impinging  at  a  zenith  angle,  f, .  The  pene¬ 
trating  sunlight  is  partially  scattered  by  the  intervening  atmospheric 
particles  in  its  passage  to  the  ground  surface.  The  scattered  light  is  dif¬ 
fused  throughout  the  atmosphere  by  further  multiple  scatterings.  The  dif¬ 
fused  light  and  sunlight  reaching  the  ground  cover  is  partially  reflected 
back  into  the  atmosphere,  thereby  generating  opposing  fluxes  of  diffused 
sky  light. 

The  observer  is  located  on  the  ground  at  the  origin  of  a  three-dimen¬ 
sional  cartesian  coordinate  system  (Figure  1).  The  target  is  assumed  to  be 
near  ground  level  and  the  observer-to-target  direction  is  along  the  posi¬ 
tive  x-axis  of  the  coordinate  system.  The  vertical  direction  defines  the 
positive  z-axis.  A  sky  point  is  determined  by  a  zenith  angle,  f  ,  measured 
from  the  vertical  z-axls,  and  the  azimuth  angle  a.  The  azimuth  angle  is 
measured  clockwise  about  the  positive  z-axis  from  the  observer-target  line. 
The  radiant  energy  in  a  differential  Interval  about  the  wavelength,  X, 
emitted  per  unit  solid  angle  from  the  sky  point  (f,B)  is  described  by  the 
function  Nx(f»  a  )dX.  This  term  is  a  function  of  the  sun  or  moon  position 
(rt»at)  in  the  sky,  atmospheric  conditions  (l.e.,  ground  haze,  cloud  cover, 
etc.)  and  the  reflectance  value  of  the  ground  cover. 


Figure  1.  Outdoor  lighting  coordinate  system. 


The  sky  light  reaching  the  observer  from  the  point  (f,a)  is  the  sum  of 
primary  scattered  light  and  multiple  scattered  light.  The  primary  scattered 
light  is  sunlight  which  has  been  scattered  toward  the  observer  by  all  at¬ 
mospheric  particles  along  the  direction  (f,a).  Let  IoX  denote  the  intensity 
of  a  collinated  beam  of  monochromatic  sunlight  at  the  top  of  the  atmos¬ 
phere.  The  intensity  of  sunlight  which  reaches  the  atmospheric  height  z  is 
given  by  (Figure  2):  ^ 

XX  "  IOX  exP  <  <sec  /  0xd*>» 

where  )3X  is  the  atmospheric  extinction  coefficient  and  exp  is  the  exponen¬ 
tial  functional  (e).  The  amount  of  this  sunlight  scattered  toward  the  ob¬ 
server  is  this  quantity  times  the  atmospheric  volume  scattering  coeffi¬ 
cient,  cr^x.  The  height  z  is  on  the  line  (f,a),  and  the  scattering  angle,  <p, 
between  this  line  and  the  impinging  sunlight  (f,,a,)  is  given  by: 

cos tp  -  cos  f#  •  cos  f  +  sin  ft*sin  f*cos(a-a#). 

This  primary  scattered  light  is  attenuated  by  a  factor 

exp  (  (sec  H  / 0xd£), 

i 

due  to  multiple  scattering  in  its  passage  to  the  observer.  The  total  pri¬ 
mary  radiance  reaching  the  observer  along  the  direction  (f,a)  is  given  by 
the  integral, 

Nx  *  Iox  J °<£XexP  (  (sec  (,)  J  0Xd£)  exp  (  (secf)  secfdz  (1) 

't  *t  1 

The  multiple  scattered  light,  on  the  other  hand,  is  diffused  light 

scattered  to  the  observer  from  the  atmosphere  (Figure  3) .  Let  SX  denote  the 
sum  of  the  flux  per  unit  area  of  diffused  light  flowing  vertically  in  both 
the  upward  and  downward  direction  at  the  atmospheric  height,  z.  The  amount 
of  flux  per  unit  solid  angle  is  0xSx/47r.  The  quantity  reaching  the  observer 
along  the  direction  (f,a)  is  this  amount  times  the  attenuation  factor: 

exp  (  (secf)  f  W) 

i 

The  total  radiance  due  to  multiple  scattering  which  reaches  the  observer 
from  all  volume  elements  along  the  direction  (f,a)  is  given  by  the  inte- 
8ra1’ 

Nx  -  f  **SX  exp  (secf*  f  0XdO  secfdz.  (2) 

*t  * 

Tousey  and  Hulburt  make  two  assumptions:  (1)  the  atmospheric  coeffi¬ 
cients  are  constant  independent  of  the  vertical  z-coordinate,  and  (2)  the 
distribution  of  diffuse  radiation  may  be  approximated  by  upward  and  down¬ 
ward  fluxes  of  radiation  along  the  z-axis.  Equation  (1)  for  the  primary 
radiation  reduces  to: 

Nx  -  Ift>a  °»X  e  ~fixztMC*s-«-flxZt  (3) 

0X  l-sec  •  cos  ^  " 
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where  the  scattering  angle,  <t>,  la  the  angular  difference  between  the  line 
If,  a)  and  the  direction  of  the  incident  sunlight  (fg»a#)«  The  equivalent 
vertical  atmosphere  thickness,  Z, ,  is  calculated  from  the  US  Standard  At¬ 
mosphere  Data  (29). 

Tousey  and  Hulburt  solve  equation  (2)  for  the  multiple-scattered  radi¬ 
ance  by  assuming  that  the  atmosphere  is  a  non-absorbing  medium,  and  consid¬ 
ering  the  boundary  conditions  on  the  vertical  fluxes  of  the  scattered  light 
at  the  top  and  bottom  of  the  atmosphere  (5,  10,  29).  Equation  (2)  for  the 
multiple-scattered  radiance  reduces  to: 

N T  -  IQX  C0SSg  /  CZT  +  (a--£-  *cosO  (l-e'^AZT,#c^ ) 

4*  1  +  SZT  \  T  0X 


+  b.  (e~^A  ZT  s*c^s  _  e-0\Zlt,ct  )  |  (4) 

( 1-sec •  cosf)  / 


The  equation  constants  are  given  by: 

a  -  2(l+gZT)k-(k-kT+T)  (1-rg) , 
b  -  (l+gZT)  (l-2k) , 
c  -  -2g(k-kT+T), 
g  -  (l-rg)  (l-n)0x» 
k  -  V  +  (l-v)  cos f #  , 

T  ■  e  _^XZTSeC^s« 

Equations  (3)  and  (4)  are  functions  of  the  volume-scattering  coeffi¬ 
cient,  cr^,  the  extinction  coefficient,  0^,  and  the  forward-scattering  co¬ 
efficient,  v .  These  atmospheric  coefficients  are  functions  of  the  type, 
size,  and  number  of  particles  in  the  atmosphere,  including  air  molecules 
and  aerosol  particles,  as  well  as  water  vapor.  The  volume-scattering  coef¬ 
ficient,  0"(£x»  determines  the  proportion  of  monochromatic  light  of  wave¬ 
length,  A,  which  is  scattered  by  a  unit  volume  of  material  in  a  direction 
at  an  angle,  <t>,  to  that  of  the  incident  light.  The  extinction  coefficient, 
is  the  total  proportion  of  light  scattered  by  a  unit  volume,  and  is 
defined  by  the  angular  integration  of  the  volume-scattering  coefficient  ac¬ 
cording  to: 

7T 

0X  ■  2tt  J'  a-^sln^d^. 

O 

The  forward-scattering  coefficient,  V ,  Is  that  proportion  scattered  in  the 
direction  of  the  incident  light. 
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Air  molecules  scatter  light  as  described  by  Rayleigh's  Law,  (22,  29) 
and  the  volume-scattering  coefficient  is  given  by: 

<r^y(air)  -  c#X~4  #  (l+k  *  cos2d>) ,  (5) 


where  k  accounts  for  polarization  defects  and  is  nearly  equal  to  unity,  and 
c  is  a  function  of  particle  density.  The  extinction  coefficient  is  given 
by: 


/3x(air) 


-4 


(6+2k)  as 


16ttc\~^ 
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Equation  (5)  shows  that  air  molecules  scatter  as  much  light  backwards  as 
forward  along  the  incident  light  direction,  and  the  forward-scattering  co¬ 
efficient  is  equal  to  y  m  0.5  for  a  pure  air  atmosphere. 

Now  let  us  consider  the  case  of  an  atmosphere  with  pure  air  and  light 
haze.  In  general,  the  atmospheric-extinction  coefficient  has  the  form  - 
c(VQ/\  )n ,  where  \.0  is  the  wavelength  of  green  light,  \0  -  0.53/x.  For  light 
haze,  c  ranges  from  0.2  to  0.1,  and  the  exponent  n  has  a  corresponding 
range  of  2.0  to  1.3.  The  ratio  of  the  wavelengths  of  green  light  to  that  of 
red  is  0.82,  and  that  of  green  to  blue  is  1.15.  The  mean  value  of  the  at¬ 
mospheric  extinction  coefficient  is  £=*0.02  km*1  (22,  29).  The  equivalent 
atmospheric  thickness  is  ZTa:8  kilometers  (29).  Note  that  for  these  values, 
the  exponential  function  ext  (-0ZTsec0)  may  be  approximated  by  l-/3ZTs ecfl 
for  values  of  8  equal  to  or  smaller  than  70  degrees.1 


Air  molecules  are  much  smaller  in  size  than  the  wavelengths  of  light 
waves.  Water  particles  tend  to  increase  in  size  and  number  with  increasing 
relative  humidity.  Assuming  spherical  particles  and  little  absorption, 
which  is  true  for  water  vapor,  the  scattering  coefficient  for  a  given  wave¬ 
length,  \,  changes  with  particle  size,  a,  in  accordance  with  Mie's  electro¬ 
magnetic  scattering  theory  (14,  22,  30).  As  particle  size  increases,  the 
angle  of  minimum  scattering  shifts  from  4>  equal  to  90  degrees  for 
Rayleigh's  Law,  backward  toward  <t>  equal  to  120  degrees.  Concurrently,  the 
scattering  distribution,  symmetrical  about  the  plane  defined  by  <t>  equal  to 
90  degrees,  shifts  forward  toward  the  incident  light  direction,  <t>  equal  to 
zero.  In  addition,  the  dependency  of  the  scattering  coefficient  upon  wave¬ 
length  changes.  Tables  listed  in  Middleton  (22)  suggest  that  the  scattering 
coefficient  may  be  approximated  by: 

°^(a/\)  -  C\n(l+k  cosa(f<«),  (35) 

where  k,  a,  and  f  tend  to  increase  with  increasing  values  of  a/\ .  The  ex¬ 
ponent  n  increases  from  -4  for  Rayleigh's  Law  to  a  small  positive  number, 
and  then  describes  a  decaying  oscillation  about  the  n-0  line  for  large 
values  of  a/\  (22). 

In  general,  the  atmospheric-scattering  coefficient  is  the  sum  of  all 
coefficients  for  the  different  particle  types.  A  similar  remark  applies  to 
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We  now  consider  Che  application  of  equations  (3)  and  (4)  derived  by 
Tousey  and  Hulburt  (29)  to  the  outdoor  lighting  problem  for  an  atmosphere 
with  pure  air  and  light  haze. 

The  equations  for  sky  radiance  obtained  by  summing  equations  (3)  and 
(4)  reduce  to: 

N\  (f’a)  "TfiT  {  3/*(1+co*2<W  +  rgC08*s  }  C08f>  (6) 


where  f,  f,— 70°.  How  is  this  expression  related  to  the  sky  irradiance,  H*, 
measured  on  a  white  diffusive  panel  placed  in  a  horizontal  position?  The 
sky  irradiance  is  given  by  the  equation: 


2k  k/2 


"'■/  / 


N^(f,a)  cosfsin fdfda, 


and  since  the  contribution  from  the  horizontal  sky  is  insignificant,  we  use 
equation  (6)  to  obtain: 


h‘“t^  r^7  (1+r8cosf8)- 


(7) 


Note  that  the  sky  irradiance  depends  upon  the  ground  reflectance,  r„  ,  and 
the  sun  elevation,  ft,  through  the  cosine  term  and  also  indirectly  from  the 

tO\  term- 

In  general,  the  sky  radiance  is  expressed  in  terms  of  the  sky  irradi¬ 
ance  by  summing  equations  (3)  and  (4)  following  reduction  for  the  atmos¬ 
phere  conditions,  and  substituting  equation  (7).  The  result  is: 


the  atmospheric  extinction  coefficient.  This  is  reflected  in  the  mean  value 
of  the  atmospheric  extinction  coefficient,  which  Increases  from  0.017  km-1 
for  pure  air  to  0.1  for  light  haze,  1  for  dense  haze  and  10  to  100  for  fog. 
The  dependency  on  wavelength  shows  a  corresponding  decrease,  becoming  inde¬ 
pendent  of  wavelength  in  heavy  fog  (22).  The  forward-scattering  coeffi¬ 
cient,  v,  shows  a  similar  shift  in  value,  ranging  from  7  -  0.5  for  pure  air 
(29,  15)  to  V  ■  0.8  for  an  extremely  turbid  medium  (21). 
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NX(C,*)  .  K  OW1 

2jr(l+rgcos£a) 


P\2Ta 


+  (l+cosC8  +  a(l-cos£)  )  (l-e"P\zTsec^)  (8) 

+  (3/4)(l+cos2<£)secCs  -  coaCa)(e  T  ^ 

1-secCgCOsC 


where 


a  -  l/2(l-rg)(l+cosCa  +  (l-coaCa)  •  e  -P\zT8ec^s)» 


The  results  lead  to  equation  (9)  when  the  exponential  functions  are 
replaced  by  their  linear  equivalents,  valid  when  both  the  sun  elevation  and 
sky  elevation  angles  are  less  than  70  degrees. 


Nx(C,a)  -  - ^ - 

2x(l+rgcosCB) 


3/4(l+cos2d>) 


+  rgcosC8 


(9) 


A  similar  expression  may  be  derived  from  equation  (8)  for  the  horizontal 
sky  radiance,  Nx  (a)  -  Nx(y,a)  at  sun  angles  C*  £70°,  by  setting  £«7r/2  and 

using  e"^ZT8ec^8  ~  1  -0xZTsecr,a. 

Nx(a)  ■  H%  ^P\ZT^  1|3/4(l+sln2^.  •  cos2(a-a#)  )+rgcos^t 
2x(l+rgcosCa) 

+  1/2PXZT  •  | l-rg  +  (l+rg)co8?a  (1°) 

-(3/2)secCg  •  (l+sin25a  •  cos2(a-ag  ))|  j 
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1 


....  > 


See  Footnote  2  for  derivation  of  the  eky  luminance  and  illuminance 
from  equations  (9)  and  (10) .  See  Footnote  3  for  derivation  of  sky  lighting 
equations  for  overcast  cloud  conditions. 


2 

The  sky  luminance,  B(£,a)  is  determined  by  the  sky  radiance  in  accordance 
with, 


B«,a) 


K  •  Nx  (? ,a)d\  , 


(36) 


where  i|»x  is  the  luminosity  function  for  the  light-adapted  eye.  Using  equa¬ 
tion  (9),  the  sky  luminance  for  £,€s  £70®,  is  given  by: 

Bf£,a)  - _ Ei _  { 3/4(l+cos2^)  +  r»cos^8lsec?  (37) 

2:r(l+rgCOS?fl)  1  ’ 


where  E,  is  the  sky  luminance  and  the  ground-reflectance  values  are  con¬ 
stant  independent  of  spectra  (i.e.,  snow  or  black  earth).  Using  equation 
(10)  in  equation  (35),  the  horizontal  sky  luminance,  Bh(a)  for  <70®  is 
given  by: 

Bh(a)  -  E«(IW _  (3/4(l+sin25  •cos2(a-a,))  +  racosC8 

2T(l+rgcos?a)  *  8 

+(1/2)PkZt  •  |  l-rg+(l+rg)cos?8.  (38) 

-(3/2)sec£g  •  (l+sin2?a  •  cos2(o-ag))|  | 

3 We  consider  for  comparison  an  atmosphere  with  a  heavily  overcast  sky.  The 
scattering  coefficient  is  dominated  by  that  for  the  water  particles  within 
the  cloud  layer.  These  droplets  are  large  in  size  (20  microns  average  dia¬ 
meter)  and  numerous  (5  to  50xl0'6  drops  per  in3)  (22,  24).  The  correspond¬ 
ing  extinction  coefficient  is  P\^10  to  100  km"*,  and  it  is  relatively  in¬ 
dependent  of  wavelength  (22).  The  scattering  process  is  dominated  by  the 
cloud  layer,  and  the  effective  atmospheric  thickness  of  the  turbid  medium, 
t.  Assuming  that  t's  magnitude  is  on  the  order  of  one  kilometer  or  larger, 
then  all  exponential  functions  tend  to  zero.  The  sky-radiance  is  determined 
by  equation  (4)  for  the  multiple-scattering  process,  and  this  equation  re¬ 
duces  to: 

Nx(C,o)  -  k  (l+rg)  {  l+2(-^4)  •  (l-H)cosC  }  (39) 

4*(l+gt)  8  1  1+rg  1 


Equation  (39)  is  a  function  of  the  zenith  angle,  £,  and  ground  reflec¬ 
tance,  r8.  The  maximum  value  occurs  at  the  sky-point  directly  overhead,  and 
the  sky-radiance  decreases  with  the  increasing  values  of  zenith  angle.  The 
equation  is  independent  of  azimuth  angle,  a,  end  assumes  a  constant  value 
independent  of  zenith  angle  when  a  ground  surface  is  highly  reflective  as 
for  a  snow  field  (rg*l).  The  sky  lrradiance  is  given  by: 
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Ground  Rad lance 


The  ground  surface  Is  assumed  to  be  covered  with  a  material  of  uniform 
composition  having  a  diffusive  reflectance,  rQ,  which  is  spectrally  depend¬ 
ent.  The  radiance,  Wx,  emitted  in  any  direction  is: 

W\  "  7  rg  H“  •  <U> 


where  H*  is  the  sky  irradiance.  The  observer-target  distance  is  assumed  to 
be  short  enough  that  the  sky-radiance  distribution  remains  the  same  over 
all  ground  points  in  keeping  with  previous  assumptions.  Consider  an  incre¬ 
mental  area  of  sky  located  along  the  line,  (£,a)  as  measured  from  a  coordi¬ 
nate  system  transposed  to  a  ground  surface  element  (Figure  4).  The  sky 
point  emits  radiance  Nx(£,«)  in  the  direction  of  the  surface  element.  Its 
angular  area  is  dw  -  sin(jd£da.  The  light  from  the  sky  point  is  incident  at 
an  angle  £  to  the  ground  surface  normal.  The  Irradiance  at  the  ground 
element  from  the  sky  point  is  dHx  (£,a)  ■  Nx(£ ,a)*cos^«sln^d^da.  The  total 
irradance  from  the  sky  is: 

f/a 


2ir  70 

H*  -  J  j  Nx  (5,a)cos^*sin?dCa 


H°  -  ha.  c°?jg.  k  (l+rR)  |l+4/3(il!i)  •  (1-7?)} 

K  1+gt  1  l+r„  ‘ 


(40) 


and  is  a  function  of  the  zenith  angle,  and  ground  reflectance,  rQ.  Sub¬ 
stituting  equation  (40)  into  (39)  results  in: 


MS.a) 


H* 

l+(4/3)f 


(l+2f  cos?). 


where  f 


1-r 


1+r 


i)  (l-  *7) 
8 


(41) 


The  sky  radiance  may  be  expressed  in  terms  of  the  horizontal  sky  radi¬ 
ance  ■  Ny(7r/2),  as  follows: 

Nx(0  -  N£(l+2f.cosD,  (42) 


where  the  parameter  f  is  the  same  as  in  equation  (41).  Assuming  that  the 
ground  reflectance  is  independent  of  spectra,  and  introducing  the  luminos¬ 
ity  function  for  the  light-adapted  eye,  the  corresponding  sky-luminance  is 
given  by: 


B(0  -  Bh(l+2f.cos0, 


(43) 


where  Bh  is  the  horizontal  sky-luminance.  At  least  in  form,  this  is  in 
agreement  with  the  experimental  measurements  by  Noon  and  Spencer  (23)  of 
sky  luminance  on  overcast  days. 


Rather  chan  solve  this  equation,  we  note  that  H*  may  be  measured  directly 
with  a  spectrometer  placed  In  a  horliontal  plana. 


Target  Radiance 

The  target  la  assumed  to  be  a  vertical,  flat-panel  emplaced  at  height, 
h, ,  above  the  ground  level  and  of  such  small  slse  that  It  does  not  cast  a 
significant  shadow  upon  the  ground  Immediately  about  It.  The  target  Is  ori¬ 
entated  so  that  the  surface  normal  is  parallel  to  the  x-axis. 

The  target  surface  Is  coated  with  a  diffusively  reflecting  material 
having  a  spectrally  dependent  reflectance  value,  rt.  The  radiance  emitted 
In  any  direction  from  the  target  surface  Is: 


T  1  T 
W\  »  x  rt  HX  • 


(12) 


where  Hx,  the  total  irradiance  at  the  target  surface,  is  the  sum  of  the 
irradiance  from  the  sky  and  that  from  the  ground  surface  to  the  target 
front;  l.a.,  h(  -  HTXS  +  h\b  . 

We  now  derive  each  of  these  lrradiances  In  turn  and  consider  their 
forms  for  the  case  of  a  pure  air  atmosphere  mixed  with  slight  haze.  The  ob¬ 
server-target  distance  Is  assumed  to  be  short  enough  that  the  sky  radiance 
distribution  remains  the  same  at  all  ground  positions  In  keeping  with  pre¬ 
vious  assumptions.  A  radiance  Nx(£,m)  le  emitted  from  the  sky  point  (£,a) 
as  measured  by  a  coordinate  system  transposed  to  the  target  position  (Fig¬ 
ure  5).  The  sky  point  has  a  differential  solid  angular  area,  dw  "  sin£d^da. 
The  light  from  the  sky  point  Is  at  an  angle  ^  to  the  surface  normal  of  the 
target,  where  coe^-  sln^«cos(x-a).  The  total  Irradiance  at  the  target  sur¬ 
face  from  the  sky  is: 

iir/l  */i 

/ 

*/l  o 

The  sky  radiance  for  pure  alr/sllght  haze  atmosphere  Is  described  by 
equation  (8).  However,  this  equation  contains  exponential  functions  which 
produce  singularities  following  integration.  For  this  reason,  we  approxi¬ 
mate  equation  (8)  by  equation  (9)  for  sun  and  sky  zenith  angles  less  than 
70  degrees;  i.e.,  £.,££70*,  and  by  equation  (10)  for  sky  angles  near  the 
horizon;  i.e.,  70*  90*.  The  target  is  a  vertical  surface  and  as  we 
shall  see,  only  the  horizontal  sky  provides  a  significant  contribution  to 
the  sky  Irradiance.  We  separate  equation  (13)  for  the  sky  irradiance  Into 
the  sum  of  two  components:  (1)  an  Integration  of  equation  (9)  over  the 
range  0<£<  70*,  and  (2)  an  Integration  of  equation  (10)  over  the  range  70* 
<£<  90*. 


Nx(£,a)  sin^tcosad^da 


(13) 


16 


V 


VT 


. . 


The  results  of  the  Integration  of  the  first  component  are: 


hTsO)  - 


—  i  0.135-(0.372raecosf8 

'f8>  ' 


2x(l+rgecos4 
+0.79cosf8*slnf8acos  f8+0.645sln2f8 
+0.231sln2f8«cos2fl8)| 

while  those  for  the  second  component,  using  ~  0*16,  are: 

hTs(2)  -  — - - - |l.59  {  3+4r0cosf.+  sin2f,*(l+co82aJ  l 

v  2ir(l+rgcosfs)  l  '  ^  S 

-0.255  |  2(l-rg  +  (l+rg)  •  cosf8  -  1.5secr8) 
f88in2  fg  •  (l+cos2a8)|  | 


-sec 


(14) 


(15) 


A  comparison  of  the  two  components  shows  that  the  sky  lrradlance  Is  domina¬ 
ted  by  equation  (15)  for  the  horizontal  sky. 


We  now  consider  the  lrradlance  received  by  the  target  from  the  ground 
(Figure  6).  The  radiance  emitted  from  an  Incremental  element  of  ground  area 
located  at  a  radius  p  and  angle  a  measured  from  the  target's  front.  Is 
W*(p,a)  given  In  equation  (11).  The  straight  line  distance  between  the 
ground  element  and  target  Is  /-(p^h2)1^  ,  and  the  ground  radiance  Is  at¬ 

tenuated  In  Its  passage  through  the  atmosphere  by  an  exponential  function 
of  the  spectral  dependent,  atmospheric  extinction  coefficient,  Further¬ 
more,  a  certain  portion  of  the  horizontal  sky  light,  N^(x/2,a),  Is  added  to 
the  ground  radiance.  The  net  radiance  at  the  target  from  the  direction  at 
the  ground  element  is  given  by: 

wf(p,a)  -  W?  (p,a)e'CT\^  +*v(>r/2  ,a)(l-e-<rV*  ).  (16) 

The  angular  area  subtended  by  the  ground  element  Is  dip  »  ds*cosd/j22, 
where  ds  *  pdadp  and  the  angle  9  Is  measured  from  the  vertical  z-axls  to 
the  straight  line  along  the  distance,  1;  l.e.,  cos®-  h/^.  The  light  from 
the  ground  element  Is  at  an  angle  ^  to  the  surface  normal  of  the  target, 
where  cos  ^  -  sln0*cos(x-a) .  The  total  lrradlance  at  the  target  from  the 
ground  la: 

in/ 2  „ 

H^g  -  -h,/ /  |l!^'N)Ua))e'a^+N*(a)^  >^2  cosadida.  (17) 

V*  ht 


Equation  (17)  contains  exponential  functions  which  produce  singulari¬ 
ties  when  Integrated.  However,  the  values  of  the  extinction  coefficient, 
crv,  are  mail  In  a  pure  air  atmosphere,  and  the  exponential  function  de¬ 
creases  at  a  slower  rate  with  Increasing  distance,  than  does  the  Inverse 
square  function,  MS.1.  In  other  words,  the  contribution  of  distinct  ground 
elements  becomes  Insignificant  with  Increasing  distance,  before  appreciable 
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atmospheric  attenuation  can  occur.  The  horizontal  sky  radiance  terms  cancel 
in  the  equation,  and  the  resulting  integration  of  the  ground  radiance  is: 


rg* 


(18) 


Background  Radiance 


The  background  visually  surround!:.*  the  target  has  a  diffusive,  spec¬ 
trally  dependent  reflectance,  rb.  It  slopes  upward  away  from  the  observer 
with  the  surface  normal  at  an  angle  J  -  V  to  the  vertical  z-axis  (Figure 
7).  The  background  incline  is  slight  ( V  nearly  equal  to  tt/2),  and  little 
shadow  is  cast  upon  the  ground  to  its  immediate  front.  It  does  not  provide 
significant  radiance  to  ground  elements.  The  observer-target  line  inter¬ 
cepts  the  background  surface  at  a  height,  hb,  above  ground  level.  The  radi¬ 
ance  emitted  from  the  background  in  the  direction  of  the  observer  is: 

(19) 

where  the  lrradiance  on  the  background,  Hx,  is  the  sum  of  the  irradiance 
from  the  sky  and  that  from  the  ground;  i.e.,  ■  Hx*  +  Hx°  .  We  again  de¬ 
rive  each  of  these  irradiances  in  turn  and  consider  their  forms  for  the 
case  of  a  pure  air/slight  haze  atmosphere. 


The  radiance,  Nx  (f,a),  is  at  an  angle  4/%  to  the  surface  normal  of  the 
background,  where  cos^,-  co8f«8iny-8inf*cosa»cosY •  The  lrradiance  from  a 
differential  element  of  angular  area  located  at  the  sky  point  is  dHx*(f,a)“ 
N^( f ,a)»sinf«C08^#d  fda.  The  total  irradiance  at  the  background  surface  from 
the  sky  is: 


// 


dH 


bs 


(f,a)  + 


JL 

2 


TT 

n 


cos~l(cosYcosa) 


dH\* 


(f i a)  (20) 


3ir 

~r 


However,  equation  (20)  was  derived  assuming  that  angle  Y  is  nearly  equal  to 
tt/2;  that  is,  that  the  surface  incline  is  slight.  For  this  reason,  we  ap¬ 
proximate  the  sky  irradiance  at  the  sloped  background  by  that  at  a  horizon¬ 
tal  surface;  i.e., 

«\8  -  HSX.  (21) 

The  radiance  from  an  incremental  ground  element  (Figure  8)  on  a  radius 
P  and  angle  a  ,  measured  from  the  surface  front,  is  Wx°(p,a)  given  in  equa¬ 
tion  (11).  The  radiance  is  at  an  angle  ^B  to  the  surface  normal,  where 
cos^  sinYcosfl-cosYsinflcosa,  and  the  angle  0  is  between  the  vertical  z-axis 
and  the  straight  line  distance, &  -( p2  +  hb2 )  1 .  The  total  ground  irradiance 

l8:  4jL  « 

H\9  ‘  hb  f  f  {(iHjrB-N£(a)).e-°Ytf  +  N*(a)}  .  (hbsinY-(i2-hb2 

?  hj)(l-ftan2Ysec2a)  ^ 

•cosY  cos 
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a)^j  da.  (22) 
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Background  lrradiance  from  the  sky. 


The  background  surface  Is  assumed  to  parallel  the  y-axls;  that  Is,  the  line 
of  steepest  ascent  Is  In  the  x-z  plane.  Note  that  the  lower  limit  of  inte¬ 
gration  of  the  variable  j2  varies  with  the  angle,  a. 


We  assume  in  equation  (22)  that:  (1)  the  exponential  functions  may  be 
ignored  due  to  moderate  extinction  in  a  pure  air  atmosphere,  and  (2)  the 
lower  limit  of  the  distance  variable,  SL ,  is  hb«tanY*seca,  since  the  angle  Y 
is  nearly  equal  to  u/2  and  the  term  tan*Yaec2a  is  much  larger  than  unity. 
Equation  (22)  reduces  to: 


\  H*r„ cosYcotY.  (23) 

Note  that  H\°  «  0  when  Y  “  it/2  as  would  be  expected  for  a  horizontal  sur¬ 
face,  but  that  H^8  approaches  infinity  as  the  angle  Y  approaches  zero  be¬ 
cause  of  the  second  assumption  above.  The  total  irradiance  at  the  back¬ 
ground  surface  is  the  sum  of  equations  (22)  and  (23). 

Hk  -  Hg(l+ -J-rgCosYcotY)  •  (24) 


Telescopic  Coating  Irradiance 

The  exterior  surface  of  the  objective  lens  of  the  observer's  telescope 
is  partially  coated  with  a  material  having  spectral  dependent  light  trans¬ 
mittance  and  diffusive  scattering  coefficients.  The  amount  of  light  scat¬ 
tered  by  the  coating  is  a  function  of  the  irradiance  at  the  exterior  sur¬ 
face  of  the  telescope.  The  determination  of  this  irradiance  is  equivalent 
in  formulation  to  that  for  the  target,  except  that  the  telescopic  front 
faces  away  from  the  observer,  instead  of  toward  him.  The  irradiance  is  the 
sum  of  that  from  the  sky  and  that  from  the  ground;  l.e.,  -  H“  +  H\°  . 

The  irradiance  from  the  sky  is: 

TT 

T 

N^(f»a)8in2f  •  cosadfda,  (25) 

while  that  from  the  ground  is: 

C 

/ 


where  hc  denotes  the  scope  height  above  ground  level. 

Equation  (25)  for  the  sky  irradiance  is  separable  into  the  sum  of  two 
components.  The  first  component  is  the  integration  of  equation  (9)  between 
the  limits  of  0°  and  70°  for  the  zenith  angle,  while  the  second  component 
is  that  of  equation  (10)  between  70°  and  90°.  The  first  component  is: 


-  N^(a))e  a^+  N^(a)j  cosadjtfda,  (26) 


..  .-y 
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H*  / 

H?8(l)  -  _ ^0.135  -  (0.37^00*^-0.79008^. 

*  2u(l+rgco8r8)'  8 

sin(8cosas  +  0.645sin2f8  +  0.231sin2ft  cos^a.,)  j,  (27) 

and  using  P^Z-j;— 0.16,  the  second  component  is: 

H*  / 

Hf 8(2)  - - ^ _  <  1.59(3+4r0cosf  +2sin2f  (l+cos2a.)) 

V  2>r(l+rgcosra)  *  '  ‘  '  * 

+0.255(2(l-rQ+(l+rg)cosf,  -3/2secr,  )  -sec f ,  •  sin2f,  (l+cos2a, ))  }  .  (28) 


Equation  (27)  for  the  telescopic  surface  is  identical  to  equation  (14) 
for  the  target  since  they  are  the  irradiance  from  the  zenith  sky.  Equations 
(15)  and  (28)  give  the  irradiances  from  the  corresponding  horizontal  skies, 
and  are  identical  in  form  except  for  the  sign  of  the  second  bracketed  term. 


The  irradiance  from  the  ground  is  given  by  equation  (26)  assuming  that 
the  exponential  functions  may  be  ignored  in  a  pure  air/slight  haze  atmos¬ 
phere.  The  result  is: 


HCXg 


_i.  u8r 

2  Hk  g 


(29) 


Equation  (29)  is  identical  to  equation  (18)  for  the  ground  Irradiance  at 
the  target,  as  would  be  expected. 


Summary  of  Radiance  Equations 

The  irradiance  and  radiance  equations  derived  for  a  pure  air/slight 
haze  atmosphere  are  listed  in  Table  1  as  a  function  of  the  emitting  sur¬ 
face. 


TABLE  1 


Radiance  Equations 


Surface 

Irradiance 

Radiance 

Ground 

—  _ 

(ID 

Target 

(14)+(15)+(18) 

(12) 

Background 

(24) 

(19) 

Telescope 

(27)+(28)+(29) 

— 

* 
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Target-Background  Radiance 


The  radiance  emitted  from  target  and  background  must  be  adjusted  for 
atmospheric  attenuation  and  addition  of  air  light,  to  determine  those 
values  which  reach  the  observer  (22,  3,  4,  6,  7,  8,  9,  11,  and  13)  (Figure 
9).  The  radiance  at  the  observer  from  the  target  and  the  Intervening  atmos¬ 
phere  is: 

wj'-  V^e'^T+N^  ( 0 )  .  (l-e"°^T) ,  (30) 

where  the  radiance  emitted  from  the  target,  W^,  is  given  by  equation  (12), 
and  the  horizontal  sky  radiance,  Ny(0),  the  observational  direction  is 
given  by  equation  (10)  with  the  azimuth  angle,  a  •  0.  The  straight  line 
distance,  /y ,  between  the  observer  and  target  is  in  the  same  units  as  the 
inverse  of  the  atmospheric  extinction  coefficient,  a^,  along  the  path  of 
observation. 


The  radiance  at  the  observer  from  the  background  is: 

W^«  •  e_<T>^b  +  N^(0)  •  (l-e"°>A),  (3!) 

where  the  background  radiance,  W^,  is  given  by  equation  (19)  and  £.b  is  the 
straight  line  distance  between  the  observer  and  the  background. 


The  observer-background  distance  (j?b)  is  a  function  of  the  observer- 
target  distance  (j£y)  and  the  terrain.  The  straight  line  connecting  the  ob¬ 
server's  telescope  and  the  target  passes  through  the  background  position. 
By  congruent  triangles, 


k 


(hb-hc) 

(  hf-hj. ) 


(32) 


where  hc,  hT,  and  hb  are  the  vertical  heights  of  the  telescope,  target  and 
background  positions  above  ground  level.  The  straight  line  distances  are 
given  by: 


k  -  (Xx  +  (hx-hc)2)1/z  , 

k  -  (Xb  +  (hb-hc)2)1/z  ,  (33) 


where  Xy  and  Xb  are  the  horizontal  observer-target  and  observer-background 
distances.  The  directional  cosines  of  the  connecting  line  are  a*  -  XT//T  & 
a,  “  (hT-hc)/j!T.  The  apparent  angular  area  of  the  target  subtended  at  the 
observer' 8  position  is, 

<4  -  ATax//z,  (34) 


where  Ax  Is  the  target-area. 
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Geographic  Coordinates 

The  observer's  coordinate  system  specified  by  the  observer-target  line 
may  be  converted  to  the  geographic  coordinates  of  the  global  world  given  by 
(a)  geographical  latitude  of  the  observer's  position,  (b)  bearing  of  the 
observer-target  line  from  true  North,  and  (c)  declination,  meridian,  and 
azimuth  of  the  sun.  Burge  (1)  lists  equations  converting  the  zenith  and 
azimuth  angles  of  the  sun  as  measured  in  the  observer's  coordinate  system 
to  appropriate  coordinates  in  the  geographic  system. 

The  zenith  angle  of  the  sun,  is  converted  by  cos  fg  *  sln/xeslnd  + 
cost*coBM»cos5,  where  the  angle  n  Is  the  geographical  latitude  of  the  ob¬ 
server's  position.  The  angle,  5,  is  the  declination  of  the  sun,  that  is, 
the  angular  distance  of  the  sun  above  or  below  the  equator  varying  with  the 
season.  The  angle,  t,  is  the  meridian  angle  measured  on  the  equator  from 
the  North  Point  to  the  circle  passing  through  the  sun.  The  meridian  angle 
may  be  computed  from  the  local  apparent  time  t, ,  by  t  »  t,  -12  hours. 

The  sun's  azimuth  angle,  a,,  is  converted  by  a,-  a0  -aNf  where  a0  is 
the  azimuth  of  the  sun  measured  from  true  North  as  given  by  sina0  -  (sint 
•  cosf)/sinft.  The  angle  aN  is  the  bearing  of  the  observer-target  line 
measured  from  true  North. 


CONCLUSION 

Equations  are  derived  for  the  sky  radiance  of  a  clear  sky  with  slight 
haze  as  a  function  of  the  sky  irradiance,  the  sun  angles  and  the  spectro- 
reflectance  of  the  ground  cover.  These  expressions  are  used  to  derive  equa¬ 
tions  for  the  radiances  emitted  from  a  panel  target  and  its  background  as  a 
function  of  their  spectro-ref lectance  values.  The  apparent  radiances  reach¬ 
ing  the  telescope  from  the  target  and  background  are  derived  by  considering 
the  intervening  atmosphere  and  the  distances  of  the  target  and  background 
from  the  observer.  The  irradiance  arriving  at  the  telescope-front  from  the 
surrounding  sky  and  ground  are  derived  to  complete  the  lighting  scenario. 
Finally,  the  observer-target  coordinate  system  is  related  to  the  geographic 
coordinate  system  of  the  global  world. 
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F40I  BUIK-MOr 


Spectral  Data  for  Daylight,  Natural  Terralna  and  Atmospheric  Attenuation 


The  apectral  data  for  daylight  lrradiance,  epectroref lectancee  of 
natural  terrain  features  and  the  atmospheric  attenuation  coefficient  a  are 
briefly  reviewed  In  this  section. 


DAYLIGHT  IRRADIANCE  SPECTRUM 

The  relative  spectral  power  from  the  total  sky  and  sun  Irradiating  a 
horizontal  surface  changes  only  slightly  over  the  visual  spectrum  during 
daylight.  This  Is  true  for  the  sun  at  altitudes  anywhere  between  two  hours 
after  sunrise  and  two  hours  before  sunset,  and  the  sky  anywhere  from  total¬ 
ly  covered  with  clouds  (overcast  sky)  to  totally  free  of  clouds  (blue  sky). 
The  ICE  standard  llluminant  D65  Is  based  on  numerous  spectroradlometrlc 
measurements  of  daylight  at  different  locations  throughout  the  world,  and 
best  represents  natural  daylight  (Judd  and  Wyszeckl,  1975  [3]). 


SPECTROREFLECTANCES  OF  NATURAL  TERRAIN  FEATURES 

The  spectral  reflectance  of  natural  terrains  has  been  extensively  In¬ 
vestigated  by  Krlnov  (1953)  (4).  As  reported  In  Pendorf  (1954)  (5),  Krlnov 
took  most  of  his  spectral  measurements  on  clear  days  around  noontime.  In 
the  case  of  horizontal  terrains  (i.e.,  snow  fields,  grass,  earth),  the 
spectrograph  was  placed  at  ground  level  on  a  tripod  and  directed  downward 
at  an  angle  of  45  degrees  and  an  azimuth  of  90  degrees  from  the  sun.  For 
vertical  terrains  such  as  trees,  the  spectrograph  was  directed  horizontally 
and  at  an  azimuth  of  between  135  degrees  to  225  degrees  from  the  sun. 

The  spectral  curves  for  snow  fields  are  relatively  flat.  The  values 
for  freshly  fallen  snow  drop  gradually  toward  the  red  portion  of  the  spec¬ 
trum.  The  curve  for  snow  covered  with  a  film  of  Ice  Is  flat,  Independent  of 
wavelength. 

The  spectral  curve  for  vegetative  formations  has  a  weak  maximum  In  the 
yellow-green  region,  and  an  upward  slope  in  the  red  portion.  The  curve  for 
coniferous  forests  In  the  winter  period  is  relatively  flat  with  a  slight 
rise  In  the  near-infrared  region.  The  curve  for  coniferous  forests  In  the 
summer  period,  dry  meadows,  and  grass  (excluding  lush  grass),  is  higher  In 
value  and  shows  a  similar  rise  In  the  near  Infrared.  The  curve  for  decidu¬ 
ous  forests  In  the  summer  period  and  all  lush  grass  has  a  weak  maximum  in 
the  yellow-green  region,  and  a  strong  rise  in  the  infrared  region.  Ripe 
field  crops  and  forests  in  autumn  show  an  upward  slope  in  the  entire  green- 
to-red  portion,  and  a  steep  slope  in  the  near  Infrared  portion. 

Measurements  from  the  ground  and  from  the  air  are  In  agreement  for 
bare  ground  areas,  such  as  soli  and  shallow  or  dense  vegetative  founda¬ 
tions.  However,  measurements  of  forests  from  the  air  are  only  one-half  to 


one-quarter  of  chOM  for  a  ground  observer  dua  to  cha  larga  shadow  araaa 
between  Individual  traaa.  For  aueh  vartieal  vagataciva  foraatlone,  tha 
ground  spectral  raflactanca  valuaa  uaad  to  ealculata  eky-radlance  should  bo 
thoaa  valuaa  aeaaurad  froa  tha  air. 


ATMOSPHERIC  ATTENUATION  COEFFICIENTS 

Tha  spactrua  variations  of  tha  ataoapharic  attenuation  coefficient 
uaad  in  this  study  are  a  aiaplif lcatlon  of  actual  data  (Barnes,  1968  [1]), 
since  tha  line  structures  of  ataoapharic  abaorptlon  due  to  water  vapor, 
carbon  dioxide,  etc.,  have  been  Ignored. 
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Computer  Program 


A  computer  program  was  written  In  Fortran  IV  language  to  compute  the 
apparent  radiances  from  the  target  and  the  background.  The  user  must 
specify  the  angular  position  of  the  sun  (moon)  relative  to  the  observer  and 
target,  the  position  of  the  target,  and  the  spectroref lectance  values  of 
the  ground  cover,  target,  and  viewed  background. 

The  program  first  calculates  the  sky  illuminance  (foot-candles),  and 
the  luminance  ( foot-lamberts)  from  the  viewed  background  at  the  surface, 
the  luminous  reflectance  value  of  the  background,  and  the  horizon  lumi¬ 
nance.  These  values  are  printed  out  for  reference.  The  program  then  com¬ 
putes  the  apparent  radiances  from  the  background  and  target. 

The  computer  program  is  attached  below,  and  a  list  of  subroutines 
follows : 


1.  TEST  -  Main  program  sets  initial  conditions  and  calls  for  com¬ 
putation  of  target  and  background  radiances. 


2.  SUNA  -  Subroutine  computes  constants  of  irradiance  and  radiance 
equations  as  a  function  of  the  sun's  elevation  and  azimuth,  and  the  back¬ 
ground  slope  angle.  Called  by  TEST. 

3.  LUMIN  -  Subroutine  computes  sky  illumination  (foot-candles), 
the  viewed  background  luminance  (foot-lamberts),  the  luminous  reflectance 
of  the  background  and  the  horizon  luminance.  These  values  are  computed  from 
the  constants  of  SUNA  and  the  spectroref lectances  of  the  ground  cover  and 
background.  Called  by  TEST. 


4.  RAPT  -  Subroutine  computes  the  radiances  emitted  by  the  target 
and  background  at  the  specified  range.  Called  by  TEST. 

5.  SKYI  -  Data  file  of  sky  Irradiance  spectrum.  Called  by  LUMIN 

and  RADT .  — 


6.  GRDS  -  Data  file  of  ground  cover  spectroref  lectance  spectrum. 
Called  by  LUMIN  and  RAPT. 

7.  TGTS  -  Data  file  of  target  spectroref lectance  spectrum.  Called 

by  RADT. 


8.  BKGD  -  Data  file  of  viewed  background  spectroref lectance  spec¬ 
trum.  Called  by  LUMIN  and  RADT . 

9.  ATTC  -  Data  file  of  horizontal  atmospheric  attenuation  coef¬ 
ficients  for  standard  clear  atmosphere,  visibility  23.3  kilometers.  Called 
by  RADT. 
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SMYTH, STHFZ,T100G, 

ACCOUNT,  HE***, 

FTNCSL,R,T) 

MAP (PART) 

LGO. 

PROGRAM  TEST (OUT PUT, TAPC3-0UTPUT) 

COMMUM/BKGD/AHGB 

COMMON /RADD/XT (79),X3(79)»SC(79> 

C  TARGET  DIAMETER  AMD  RANGE  IH  METERS 
DATA  DTO/3,/ 

DATA  R/ 2000, / 

C  VIEWED  BACKGROUND  GO-DEGRCCS  FROM  GROUND  NORMAL  (10-DEGREt  SLOPE) 

DATA  ANGB/I, 3^626/ 

C  SUN  DIRECTLY  OVERHEAD 
DATA  EL, A2/2*0, / 

WRITE (3, 990) 

090  FOP.MATdM  , 6X,  U'IS IEILITY  COMPUTATIONS*) 

C  SKY  LIGHTING  CONDITIONS 

CALL  SUHA(EL,AZ, AHGB) 

CALL  LUHIH(ES,6L,BR»BMS) 

WRITEC3,  091)  ES,  BL,  BP.,BHS 

991  FORMAT (1H  ,6X,«SKY  1UCFC) «*,F10.%. •  BKGO  LUM( FL )• •# FlO. %, •  BKGD  R 
CEFLECTANCE-',F10.%,'  HORIZON  (FL)-»,F10.%> 

C  SET  UP  TARGET  ANO  BACKGROUND  RADIANCE  FOR  SPECIFIED  RANGE 
10  RXP«R*0T0*TANUNGB);2. 

WRITE ( 3, 9o%) R, RX8 

99%  F0RMATC2X, 'TARGET-RANGE*, 2X,F10.%,2X, 'BACKGROUND  RANGE *,2X,F 10. %) 
CALL  RADT(R,RXBI 

C  WRITE  OUT  TARGET  APPARENT  RADIANCE 
WRITE(  3, 1000)  (XT(I)»I-1»79) 

1000  FORMAT (2X,10(F10«%,1X)) 

C  WRITE  OUT  BACKGROUND  RADIANCE 

WRITEC 3, 1000) (XBtl 1,1-1, 79) 

C  WRITE  OUT  IRRADIAMCC  AT  SCOPE  FRONT 
WRITEC 3, 1000)<SC(I)»I-1, 79) 

STOP 


SUBROUTINE  SUNAi EL# AZ# BN) 

C  CONFUTES  CONSTANTS  OF  SKV  RADIANCE  EQUATIONS 

COHNOM/C OUST/ A, AB# BT#CT# DT# ST# BC#  CC» OC# EC#  BH# C H>  DH#  EH 
DATA  PI/3.14159/ 

A-COS(EL) 

AB-COS(BNJ*COT{BtU/2. 

BT«.1359-.7963*C0S(EL)*SIN<EL>*C05«AZ>-.6462*(SINIEL><'*2>-.2313*(  ( 
QSIN(EL)*C0S(AZ))**21-. 2549*{  2.*<  l.*CO$(EL>-1.5*SEC<EL)  )-SECCEL)*( 
Q5IN<EL)**2)*(l.»(COS(AZ))**2)) 

CT-.25AO*(3.  +  2.*(SIM<EL}**2)*a.*{COS(AZ)**2»n 
DT— (,372*C0S<EL)-.5*(  I.-CQE  tELM ) 

ET-A 

BC.#1359*.7963*C0S<CL)*SIM(EL>*C0S(AZ)-.6462*t (SINCEL) »**2)-.2313* 
D((SIN(EL)*CUS(AZ)) **2) ♦.  2549*1 2»*(1 .♦COS (tL)“1.5*S EC (EL) )~SEC ( EL ) * 
0(<SIN{EL))**2)*(l.^CaS(AZn**2) ) 

CC-CT 

DC  —  (,B72*C0S(CL)*.5*(1.K:0S(EL)  )) 

EC-ET 

BH..5*(l.+COS(£l)-1.5*SEC<CL>*(l.MSIN(EL>*CQS<AZ))**2>) 

CH«« 7? *(].«♦( SI M ( EL  )*C 05 tAZ))**2) 

PH-(A-1, J/2. 

EH-A 

WRITE(3#1000)A#AB#BT#CT#DT#ET#BC#CC#  DC#EC#BH#CH, DH#EH 
1000  FORMAT ( 2X#  * SKY  IRRADIANCE  CONSTANTS' /2(2X, ?( F10. A# 2X )/) ) 

RETURN 

END 

SUBROUTINE  LUMIIKES#  BL#BR# BHS) 

COHHaN/r.AD/SH,TC#r.T#r.G»SC#AT#RB 

COMMON/C OUST/ A#  A B#  BT#CT#  DT,  ET#BC#CC#  DC#  EC#  BH#  C  H»  DH#  £H 
DATA  PI/3.14159/ 

C  ATTENUATION  WITH  SLIGHT  HAZE#  STANDARD  ATHOPHERIC  HEIGHT 
DATA  CM# XN#ZT/. 00015,1.6, 8000*/ 

C  RADIANT  FREQUENCY  RESPONSE  RANGE#  VISUAL  FUNCTION 
DATA  WO#  DW/375.#  5. / 

W-WO 

ES-O. 

BL-O. 

BHS-O.I 

C  IRRADIANCE#  MICROWATTS/CH  SQ/5  NN 
DO  10  1*1,79 
W-WfDW 

CALL  STDBS( I#X#Y#Z)  • 

CALL  GP.DS(I#RG) 

CALL  BKGD{ I#  RB) 

CALL  SKYI(I) 

ES»ES+$H*Y*OW 

BL-BL»SH*PB*Y*DW/PI 

AK»CN*((530,/W)**XN) 

AV-l./(ZT*AKl 

BHS"BHS^SH*Y*OW* ( BH+CH*AW+ ( DH*EH*AW)  *RG ) /( 2. *PI* ( 1. ♦ A*RG ) ) 

10  CONTINUE 

C  ILLUMINANCE#  FOUTCANDLES 
ES«ES*.6317 

C  BRIGHTNESS#  FOOT  LAHBEP.TS  (MILLILAMBERT5 ) 

BH$»BHS*1.9649 
BL “BL* l. 9049 
BR  BL*PI/ES 
RETURN 
END 
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SUBROUTINE  RADT<RX#P.XB) 

C  COMPUTES  RADIANT  ENERGY  FROM  TARGET  AND  BACKGROUND 
COMMON/RAD/SH#  TC#P.T#  P.G»$C#  AT#  RB 

COMHON/CONST/A#AB#BT#CT#DT#ET#BC#CC#DC#EC»Brt#CH»DH#EH 

C0MHr)N/RADD/WT{79)#WB(79)#SS(791 

DATA  PI/3. 14159/ 

C  ATTENUATION  WITH  SLIGHT  HAZE#  STAND ARO  ATMOPHERIC  HEIGHT 
DATA  CN#  XN# ZT/ #00015 #1,6#  6QQ0. / 

C  RADIANT  FREQUENCY#  10  NM  INTERVALS#  380  MM  TO  770  NM 
DATA  H0#DW/370.#10./ 
u-wo 

DO  10  I-l#79#2 
W»W*DW 

CALL  GF.DSd# RGi 
CALL  BKGD(IiRB) 

CALL  TGTS(I#PT) 

CALL  SKYKI) 

CALL  ATTC (I# AT) 

AK-CNO((530./W)**XN) 

AW*1./ ( ZT*AK) 

C  IRP.ADIANCE 

ST-SH*fPl*RG*(BT*AW*CT+RG*(DT*AW*ET) )/(l.*RG*A)> /<2.*PI> 
SB-SH*(1.+RG*AB) 

SN»SH*(  BH«,CH*AW*RG*{  DH+AW*EH) ) / ( 2. *P I*  C l.*RG*A ) ) 
SS(I)-SH*(PI*RG^(BC+AU*CC«-RG*CDC*AW*EC>)/(1.*P.G*AI>/(2,*PI» 

C  RADIANCES 

WTT-3T*P.T/PI 
WBB-SB*RB/PI 
C  APPARANT  RADIANCES 

WT(I)-WTT*EXP(-AT*RX)*SN*(l,KEXP(-AT*P.Xn 

UB<I)-WBO*tXP(-AT*RXB)*SN*(l.-EXP(-AT*RXB>> 

10  CONTINUE 
RETURN 
END 

SURP.OUTIMC  I  KYI  (I) 

COMMON/ FAD/  SH»  TC#  RT#  P.G#  SC#  AT#  P.B 
DIMENSION  S € 793 

C  CIE  STANDARD  ILLUMINANT#  0-65 

C  TOTAL  SKY  IRRADIANCE#  HORIZONTAL  PANEL#  FROM  DAY-SKY  AMD  SUN 
C  RELATIVE  POWER 

C  MW/CM2/NM  AT  5  NM  INTERVALS#  FRON  380  NM  TO  770  NM 

DATA  S/50.»52.3# 54,6#6a.?»82,0,67.1#9i.5#92.5#93.4#9O.l»e6,7#95.8» 
A 10 4 ,9# 110.9# 117. # 117.4# 117.8# 11 6* 3# 114* 9# 115 .4# 115*  9# 112*4# 106. 8# 
8109*1# 109*4# 10B*  6# 107. 8# 106*  3# 104*  8# 106*  2# 107*  7# 106*  # 104*  4# 104. 2# 

C 104.#  102*#  100.  #  48*  2#  96.  3#  96*  1#  95*8#  92*  2#  88*  7#  89.  3#  90.#  89.  8#  69.  6# 
D88,6#6  7.l7#85«5#8  3.  3#B3.?#B3.I7#81.9#80.#80.1#80.2#81.2#82.3»P0.3# 
E78,3#74»!#69»7#70,7#71»6#  73 .#74. 3# 68* #61* 6# 65 *7# 69. 9#  72. 5# 75*1# 

F69, 3# 63. 6#? 5, #46, 4# 56. 6# 66. 8/ 

C  RELATIVE  FACTOR#  1.7X 
SH-1.7*S(U 
RETURN 
END 
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SUBROUTINE  GRDS ( I#  RG) 

C  SPECTRAL  REFLECTANCE  VALUES  FOR  GROUND  COVER 
DIMENSION  R(79J 

C  SOURCE  PENNDORF  1956  VEGETATIVE  TYPE  C-1D  FOREST  AUTUMN 

DATA  R/.047# .048,. 049# ,05» .051, .052* .053# .05 A# .056# .058# .06# 

A 36.061# .062# .063# • 06 A, .067# .069# .062# .076# .076# .079# .081# .083# 
6.O87#«09 1#. 097# .112# . 123#. 13 A# .142#* IF 1# .159#. 168# .173# .176# .ISA# 

C.  191# .0.92# .193# .19 5#  3*.  196#  .193#  .191#  .192#  3*.  193# .192#  .19# .191# 

D.  193#  ,20  2#  .21 2#  .223# .234# .246# *2 59# .272# .265# » 3# .315# .33# .345# 

E. 361# .378# •  393#  •  409#  . 426#  . 443#  .4  52#,  46#  .  47#  .  481#  •  491#  .  5/ 

RG-R(I) 

RETURN 

END 

SUBROUTINE  TGTS(I#RT) 

DIMENSION  P  ( 79) 

C  TARGET  SPECTRAL  REFLECTANCE  VALUES#  LUSTERLESS  O.D.  PAINT  34087 

DATA  R/5*,0409# .04 14#. 0427# .0436# 2*. 0445#. 04 5# .0454# .0464#. 0473# 
1.0479# .0491# .0509# .0527# .0536# .0554# .0564# .06# .0618# . 0636# .0645# 
2.0682#  .072% ,0745#. 077 3#  .08#. 0609#  .0818#  46.0627#  46.0836#  26.0827# 

3. 08  36# .0854# .0891#. 0909# • 0927# . 0945# . 0973# . 0991# , 1# . 1009# .1018# 
4.1054# .0064#. 106 2# .1091# .11# .1145# .1173# .119#. 1236# .1273#. 1318# 
5.1364# .1436, .1473# .1527# .1554#. 1636, .1673# .1736# . 1816# .1909, .2# 
6.2091#  .2182#  .  2273#  .2364/ 

RT-RCI) 

RETURN 

ENO 

SUBROUTINE  BKGD( I#  RB ) 

C  SPECTRAL  REFLECTANCE  VALUES  FOR  TARGET  BACKGROUND 
DIMENSION  R ( 79) 

C  SOURCE  KRINOV  1953  TYPE  III  FIELD  GRASS  AUTUHN 

DATA  R/.03#, 0305#. 031#. 0315# ,032#.0325#a 033# .034 5# .036# .0375# .039# 
U. 0405, .042# .0425# • 043# .044# ,045# .046# .047# .048#. 049# .0495# .05# .051 
U# ,052# .054# , 056# . 059# . 062, , 0665# . 071# . 074# . 077#. 0785# • 08# . 0805# 

U. 081# .0805# ,06# .079# .076 #.077#. 076# .0765# .077#, 078# .079# .079 5# ,08# 
U. 061, .08 2# .082# .082# .081 #.08#. 081# .082# .083# • 084#. 092# .099#. 106# 

U.  116#  .026#  .14#  ,154#  ,  167#  .178#  .189#  .203#.  216#  .236#.  259#  .267#  .27  5# 
U.285#.295#«  323# , 35/ 

RB-R(I) 

RETURN 

END 


SUBROUTINE  ATTC( I, AT) 

C  ATMOSPHERIC  HORIZONTAL  ATTENUATION  COEFFICIENT  (1/K.H) 

DIMENSION  A<79) 

C  VISIBILITY  23.5  KM  STANOARD  CLEAR 

DATA  A/. 254,  .2505,  .247,.2435,,24,.  2365/  .233#  .2295,, 226,.  2225#  .219, 

1.2155. . 212/. 2065/. 205#. 203,. 201,. 199,. 197#. 195/. 193,. 191,. 189,. 1*7 

2..  185. .183. .101. .179. .177. .175. .173. .171. .169. .167. .165. .1636. .162 

36..  161  A,  .1602,  .‘159,  ,  15  7g,  ,1566, .1554, .1542,  .15  3,, 15235, .1517, .1510 
<•5,  .1504,.  14975,  . 1491, . 14845, . 1478, , 14715, . 1465,. 1458 5, . 1452, . 14455 

5,  .'1439  0,  •  1A325,  .  1A26,  •  1A195,  .  1A13,  •  1A065,  •  1A,  ,  13 95, ,139,  .138 5,  .136 

6,  .  1375,,  !37,,1365,.  136,.  (1355,,  135,.  1345,.  13A,.  1335,.  133/ 
AT-AO/IOOO. 

RETURN 

END 

SUBROUTINE  STUBS < I, XS, YS, IS) 

C  STAHOARD  OBSERVER  Cl E  1931  ONE-DEGREE  TGT 
C  380NM  TO  770HM  IN  5HH  IllCf.CHEMTS 
DIMENSION  X(79),Y(79),2(79) 

DATA  X/ .0014, .00 22, .0042, .0076, .01 43, .0232, .04 35, .07 76, .1344, .21 4 8 

1. . 2639.. 3285.. 3483.. 3481.. 3362.. 3187.. 2908.. 2511. .195A,. 1A21,. 0956 

2. . 058. .032. .01A7,.00A9,.O02A,, 0093, .0291, .0633,. 1096,. 1655, .2257, 
3. 2  9  04,. 3597,. 43  3  4,. 5  121,. 5  945,. 6  78 4,. 76  21,. 8425,. 9163,. 9786,  1.0263 
4,1. 0567, 1.0622, 1,0456# 1.0026,, 9384,. 8544,. 7514,, 6424,. 5419#. 4479, 

5. 3608..  2835,,  218*',  .16A9,  .1212,  .0874,  .0636,. 0A68,. 0329,. 0227,. 0158, 
6,0114,  .0081, .0058, .0041, .0029, .002, . 0014, . 0010,. 0007, . 0005,. 0003, 
7.0002, . 0002, . 0001, .0001/ 

DATA  Y/.0,,0001, .0001, .0002, ,0004, ,0006, .0012/.0022, .0040, .0073, 

I. 0116,. 0168,. 023,. 0298,. 038,. 048,. 060, .0739,. 091,. 1126, .139,, 1693, 
2.208,  .2586,  .  323,  .4073,  .503,  .  6082,  .  71,  .  7932,  ,  862,  ,9149,  .  954,  .  9803, 
3.995,1. 0002,, 99S, .9786,. 952, .9154, .87, .816 3,. 757, .6949, ,631, ,5668, 

4.503..  4412.. 381. .321. .265. .217. .175, *1382, .107,, 0816,. 061, .0446, 

5. 032,. 0232, .017, .0119, . 0082, ,0057, .0041, .0029, .0021, .0015, .001, 
6.0007,  .0005,  .0004,  .0003,  .0002,  . 0001,  .0001, •  0001,  2*0. / 

DATA  Z/.  10065,  .0135,.  0201,  •  0362,  .0679, .1102,,  2074, .3713,.  6456, 1.039 

II,  1.3856,1.623,1.7471, 1.7826,1.7721, 1.7441,1.6692,1.5281,1.2876, 
21. 0419,, 613, .6162,. 4652,. 3533,. 2^2, .2123, .15 82,. 1117, .0782,. 0573, 
3.0422,. 0298,  ,020 3,  .0134,  .  0067,  .  0057,  .0039,  .0027,  .0021,  .0018,  .0017, 
4.0014, .0011, , 001, . 0008, . 0006, , 0003, .10002, • 0002, . 0001, 29*0. / 

XS-X(I) 

YS-Y(I) 

ZS-Z(I) 

RETURN 

END 


LIST  OF  NOTATIONS  AND  SYMBOLS 


spectro-ref lectance  of  the  terrain  ground  cover 

spectro-ref lectance  of  the  target  surface 

spectro-ref lectance  of  background 

spectro-transmittance  of  the  coating  material 

spectro-scattering  coefficient  of  the  coating  material 

coodinates  of  a  right-hand  cartesian  coordinate  system  located  at  the 
observer  with  the  x-axis  in  the  target  direction  and  the  z-axis  ver¬ 
tically  toward  the  sky 

thickness  of  the  atmosphere 

zenith  angle  measured  from  the  z-axis 

azimuth  angle  measured  from  the  x-axis  in  a  clockwise  direction  about 
the  positive  z-axis 

zenith  angle  of  the  sun 

azimuth  angle  of  the  sun 

wavelength  of  monochromatic  radiant  energy 

intensity  of  sunlight  at  the  ozone  layer 

sky  radiance 

atmospheric  extinction  coefficient 

atmospheric  volume  scattering  coefficient 

scattering  angle 

forward  scattering  coefficient 

primary  sky  radiance 

multiple  scattered  sky  radiance 

sky  irradiance 

horizontal  sky  radiance 


43 


(Continued) 


PffflCKHNO  FAOS  BU1K-M0V  YILMR 


wx* 

ground 

radiance 

wx* 

target 

radiance 

Hx* 

target 

irradiance 

htxs. 

target 

irradiance 

from 

sky 

„Tg 

Hx  » 

target 

irradiance 

from 

ground 

p . 

ground 

radius 

/,  straight  line  distance 
hT,  target  height  above  ground  level  (x-y  plane) 
cr^,  ground  level  atmosphere  extinction  coefficient 
do),  differential  element  of  angular  area 
ds,  differential  element  of  ground  area 
$,  angle  between  z-axis  and  radiating  point 

angle  between  surface  normal  and  radiating  point 

n,  surface  normal 

Y,  angle  between  incline  of  background  surface  and  z-axis 

W^,  background  radiance 

H^,  background  irradiance 

background  Irradiance  from  sky 

H^8,  background  irradiance  from  ground 

hb,  height  of  background  on  observer-target  line  above  ground  level  (x-y 
plane) 

'pt,  angle  between  surface  normal  and  radiating  sky  point 

^gi  angle  between  surface  normal  and  radiating  ground  point 

H^,  coating  irradiance 

Hc^*  ,  coating  irradiance  from  sky 
Hc^  ,  coating  irradiance  from  ground 

(Continued) 
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V 

X 


b‘ 

X 


T 

b 


T 

I 

b 

« 


T; 
AT 
M 
6 
t 
t  a 


height  of  coating  above  ground  level  (x-y  plane) 

apparent  radiance  at  teleacope  from  target  through  intervening 
atmosphere 

apparent  radiance  at  telescope  from  background  through  Intervening 
atmosphere 

x-axis  position  of  target 

x-axis  position  of  interception  point  of  background  and  observer- 
target  line 

straight  line  distance  between  observer  and  target 
straight  line  distance  between  observer  and  background 
directional  cosine  along  x-axis  of  observer-target  line 
directional  cosine  along  z-axis  of  observer-target  line 
apparent  angular  area  subtended  by  target  at  observer's  position 
area  of  target 

geographical  latitude  of  observer  in  global  coordinate  system 
declination  of  the  sun 
meridian  angle 
local  apparent  time 

e 

bearing  of  observer-target  line  measured  from  true  North 
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